Multiphoton antiresonance in large-spin systems 
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We study nonlinear response of a spin S > 1/2 with easy-axis anisotropy. The response displays 
sharp dips or peaks when the modulation frequency is adiabatically swept through multiphoton 
resonance. The effect is a consequence of a special symmetry of the spin dynamics in a magnetic 
field for the anisotropy energy cx . The occurrence of the dips or peaks is determined by the spin 
state. Their shape strongly depends on the modulation amplitude. Higher-order anisotropy breaks 
the symmetry, leading to sharp steps in the response as function of frequency. The results bear on 
the dynamics of molecular magnets in a static magnetic field. 
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I. INTRODUCTION 

Large-spin systems have been attracting much atten- 
tion recently. Examples are S = 3/2 and S — 5/2 Mn im- 
purities in semiconductors and Mn- and Fe-based molec- 
ular magnets with electron spin S* = 10 and higher. Nu- 
clear spins 7 = 3/2 have been also studied, and radiation- 
induced quantum coherence between the spin levels was 
observed [ij . An important feature of large-spin systems 
is that their energy levels may be almost equidistant. A 
familiar example is spins in a strong magnetic field in the 
case of a relatively small magnetic anisotropy, where the 
interlevel distance is determined primarily by the Larmor 
frequency. Another example is low- lying levels of large- 5 
molecular magnets for small tunneling. As a consequence 
of the structure of the energy spectrum, external modu- 
lation can be close to resonance with many transitions at 
a time. This should lead to coherent nonlinear resonant 
effects that have no analog in two-level systems. 

The effects of a strong resonant field on systems with 
nearly equidistant energy levels have been studied for 
weakly nonlinear oscillators. These studies concern both 
coherent effects, which occur without dissipation 
and incoherent effects, in particular those related to the 
oscillator bistability and transitions between coexisting 
stable states of forced vibrations. In the absence of dis- 
sipation, a nonlinear oscillator may display multiphoton 
antiresonance in which the susceptibility displays a dip 
or a peak as a function of modulation frequency Q . 

In the present paper we study resonantly modulated 
spin systems with S > 1/2. Of primary interest are sys- 
tems with uniaxial magnetic anisotropy, with the leading 
term in the anisotropy energy of the form oi—DS^/2. We 
show that the coherent response of such spin systems dis- 
plays peaks or dips when the modulation frequency adi- 
abatically passes through multiphoton resonances. The 
effect is nonperturbative in the field amplitude. It is 
related to the special conformal property of the spin dy- 
namics in the semiclassical limit. It should be noted that 
the occurrence of antiresonance for a spin does not fol- 
low from the results for the oscillator. A spin can be 
mapped onto a system of two oscillators rather than one; 
the transition matrix elements for a spin and an oscillator 



are different as are also the energy spectra. 

We show that the coherent response of a spin is sen- 
sitive to terms of higher order in Sz in the anisotropy 
energy. In addition, there is a close relation between the 
problem of resonant high-frequency response of a spin 
and the problem of static spin polarization transverse to 
the easy axis. Spin dynamics in a static magnetic field 
has been extensively studied both theoretically and ex- 
perimentally % 0, i, i, M, EH, El m, Q d- One of 
the puzzling observations on magnetization switching in 
molecular magnets, which remained unexplained except 
for the low-order perturbation theory, is that the longi- 
tudinal magnetic field at which the switching occurs is 
independent of the transverse magnetic field [l3l. The 
analysis presented below provides an explanation which 
is nonperturbative in the transverse field and also pre- 
dicts the occurrence of peaks or dips in the static po- 
larization transverse to the easy axis as the longitudinal 
magnetic field is swept through resonance. 




FIG. 1: (Color online). Three-photon resonance in a 5 = 2 
system in the limit of a weak ac field, (a) Spin energy levels 
Em and n-photon energies nhuip- (b) Quasienergies in the 
limit of zero modulation amplitude, e^'^\m) = Em — mhuip; 
the pairwise degenerate levels correspond to one- and three- 
photon resonance, respectively. 

The onset of strong nonlinearity of the response due to 
near equidistance of the energy levels can be inferred from 
Fig. [ija). It presents a sketch of the Zeeman levels of a 
spin Em {—S < m < S) in a strong magnetic field along 
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the easy magnetization axis z. The spin Hamiltonian is 



Hq — LUqSz 



\DSl 



(1) 



where ujq is the Larmor frequency. For comparatively 
weak anisotropy, DS <C wo, the interlevel distances 
Em+i — Em are close to each other and change linearly 
with m. 

A transverse periodic field leads to transitions between 
neighboring levels. An interesting situation occurs if the 
field frequency is close to loq and there is multipho- 
ton resonance in the mth state: Nujf coincides with the 
energy difference Em+N — Em, N > 1. The amplitude of 
the resonant iV-photon transition in this case is compar- 
atively large, because the transition goes via N sequen- 
tial one-photon virtual transitions which are all almost 
resonant. Therefore one should expect a comparatively 
strong multiphoton Rabi splitting already for a moder- 
ately strong field. 

A far less obvious effect occurs in the coherent response 
of the system, that is in the magnetization at the modula- 
tion frequency or, equivalently, the susceptibility. As we 
show, the expectation value of the susceptibility displays 
sharp spikes at multiphoton resonance. The shape of the 
spikes very strongly depends on the field amplitude. 

The paper is organized as follows. In Sec. II we study 
the quasienergy spectrum and the response of a spin with 
quadratic in Sz anisotropy energy. We show that, at mul- 
tiphoton resonance, not only multiple quasienergy levels 
are crossing pairwise, but the susceptibilities in the res- 
onating states are also crossing. In Sec. Ill we show that 
multiphoton transitions, along with level repulsion, lead 
to the onset of spikes in the susceptibility and find the 
shape and amplitude of the spikes as functions of fre- 
quency and amplitude of the resonant field. In Sec. IV we 
present a WKB analysis of spin dynamics, which explains 
the simultaneous crossing of quasienergy levels and the 
susceptibilities beyond perturbation theory in the field 
amplitude. In Sec. V the role of terms of higher order 
in Sz in the anisotropy energy is considered. Section VI 
contains concluding remarks. 



II. LOW-FIELD SUSCEPTIBILITY CROSSING 
A. The quasienergy spectrum 

We first consider a spin with Hamiltonian Hq ([T]), 
which is additionally modulated by an almost resonant 
ac field. The modulation can be described by adding to 
Hq the term —SxAcosLUpt, where A characterizes the 
amplitude of the ac field. As mentioned above, we as- 
sume that the field frequency is close to luq and that 
ujf,ujo > D,A, \ujF - L0o\- 

It is convenient to describe the modulated system in 
the quasienergy, or Floquet representation. The Flo- 
quet eigenstates |V'e(^)) have the property \tps{t + Tp)) = 
exp(— i£Ti?)|i/'g(i)), where Tp = 2Tr/ujp is the modula- 



tion period and e is quasienergy. For resonant modu- 
lation, quasienergy states can be found by changing to 
the rotating frame using the canonical transformation 
U{t) = exp^—iLUpSzt). In the rotating wave approxi- 
mation the transformed Hamiltonian is 



H 



-SujSz 

5lO = LOp — LOq. 



hASx 



(2) 



Here we disregarded fast-oscillating terms 
cx Aexp{±.2iujpt). 

The Hamiltonian H has a familiar form of the Hamil- 
tonian of a spin in a scaled static magnetic field with 
components 5u> and A/2 along the z and x axes, respec- 
tively. Much theoretical work has been done on spin 
dynamics described by this Hamiltonian in the context 
of molecular magnets. 

The eigenvalues of H give quasienergies of the modu- 
lated spin. In the weak modulating field limit, A — > 0, 
the quasienergies are shown in Figlljb). In this limit 
spin states are the Zeeman states, i.e., the eigenstates 
|m)'^°) of Sz, with —S < m < S. The interesting feature 
of the spectrum, which is characteristic of the magnetic 
anisotropy of the form DSz, is that several states become 
simultaneously degenerate pairwise for A = [l^, [Tit . 
From Eq. the quasienergies e(°)(m) and e*^°)(m-f iV) 
are degenerate if the modulation frequency is 

Su) = Suj,n;N, Suj,n;N = -D + . (3) 

The condition ^ is simultaneously met for all pairs of 
states with given 2m + N. It coincides with the condition 
of A'^-photon resonance Em+N — Em = Nup. In what 
follows N can be positive and negative. There are 45* — 1 
frequency values that satisfy the condition ([3]) for a given 
S. 

The field oc A leads to transitions between the states 
|m)^°) and to quasienergy splitting. The level splitting 
for the Hamiltonian ^ was calculated earlier [ll|. For 
multiphoton resonance, it is equal to twice the multipho- 
ton Rabi frequency Qii{m; N), 



nR{m;N) = \A/2D\^^^ \D 



(5 + m + A)!(S'-to)! 



|-| ^sgn^' 



{S + my.{S - m - Ny. 



2(|A|-1)!2 



(4) 



The A-photon Rabi frequency (|3|) is cx yl'^l , as expected. 
We note that the amplitude A is scaled by the anisotropy 
parameter D, which characterizes the nonequidistance of 
the energy levels and is much smaller than the Larmor 
frequency. Therefore Qfi becomes comparatively large 
already for moderately weak fields A ^ D. 

We denote the true quasienergy states as with in- 
teger or half-integer v such that ~S < v < S. The 
quasienergies e^, do not cross. One can enumerate the 
states 1 1/) by thinking of them as the adiabatic states for 
slowly increasing (5w, starting from large negative Suj. For 
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—Sco/DS ^ 1, l^l/D the states {v) are very close to the 
Zeeman states | v) ^^'^ , with v being the eigenvalue of . 
This then specifies the values of v for all Suj. 

If the field is weak, the states are close to the cor- 
responding Zeeman states, \v) « |m)^°\ for all Suj except 
for narrow vicinities of the resonant values 6ujm;N given 
by Eq. The relation between the numbers v and m 
for w Im)^"^ is 

v = m+y 9[5uj — SLLjm-N)sgnN, (5) 

where N runs from —5* — m to 5* — m; the term = 
is eliminated, which is indicated by the prime over the 
sum; 9{x) is the step function. In obtaining Eq. ([5]) we 
took into account that, for weak fields, only neighboring 
quasienergy levels e^, and e,^±i come close to each other. 
Eq. [5]) defines the state enumerating function m{v). 
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FIG. 2: (Color online). Quasienergy levels Si, for a spin-2 
system as functions of detuning 5uj/D for the scaled field 
amplitude A/D = 0.3. The levels anticross pairwise at 
multiphoton resonances given by Eq. Q. The unperturbed 
quasienergies (the limit A = 0) correspond to straight lines 
£(°)(m) = -Slutti - Dm^/2. 

The enumeration scheme and the avoided crossing of 
the quasienergy levels are illustrated in Fig. (5] For the 
chosen S* = 2 the anticrossing occurs for 7 frequency 
values, as follows from Eq. ([3]). The magnitude of the 
splitting strongly depends on N: the largest splitting 
occurs for one-photon transitions. It is also obvious from 
Fig. [2] that several levels experience anticrossing for the 
same modulation frequency. 

B. Susceptibility and quasienergy crossing 

Of central interest to us it the nonlinear susceptibility 
of the spin. We define the dimensionless susceptibility Xv 
in the quasienergy state \ v) as the ratio of the expectation 
value of the appropriately scaled magnetization at the 
modulation frequency to the modulation amplitude, 

Xu{u:f) = {v\S-\^)IA. (6) 



In the weak field limit, Q. 

. , to(2(5cj + Dm) + DS'(S' + 1) 
^"^^ = A(R ^ n ^2 b ^ 

where m and v are related by Eq. ((51); in fact, Eq. ([7]) 
gives the susceptibility in the perturbed to first order in 
A Zeeman state |m)''°-'. 

A remarkable feature of Eq. ([7]) is the susceptibility 
crossing at multiphoton resonance. The susceptibilities 
in Zeeman states |m)*^°^ and \m + N)^^^ are equal where 
the unperturbed quasienergies of these states are equal, 
£(o)(j7j) = e^") (m-t-TV), i.e., where the frequency detuning 
is bio = Sujrn:N- In terms of the adiabatic states for 
such Suj we have from Eqs. (O, ([T]) Xi^i^p) = Xi^'i^p) for 
h'' = 1/ + sgnN. 

A direct calculation shows that simultaneous crossing 
of the susceptibilities and quasienergies occurs also in the 
fourth order of the perturbation theory provided > 3. 
Numerical diagonalization of the Hamiltonian ^ indi- 
cates that it persists in higher orders, too, until level re- 
pulsion due to multiphoton Rabi oscillations comes into 
play. 

The susceptibility Xi' is immediately related to the 
field dependence of the quasienergy e^. Since {i/\S+\u) = 
(i'|S'_|i^), from the explicit form of the Hamiltonian ([2]) 
we have 

Xu = ~2A-^de,/dA, (8) 

Simultaneous crossing of the susceptibilities and 
quasienergies means that, for an A^-photon resonance, 
the Stark shift of resonating states is the same up to 
order N — I in A; only in the iVth order the levels 
and eiz+sgnAf become split [by 2nR{m; N)]. Respectively, 
the susceptibilities Xi/ and Xi'-i-sgnAr coincide up to terms 
cx A'^I"'^. The physical mechanism of this special be- 
havior is related to the conformal property of the spin 
dynamics, as explained in Sec. IIVI 

Equation ([7]) does not apply in the case of one-photon 
resonance, = 1: it gives |xi/| oo for Slu (5wm;i- 
This is similar to the case of one-photon resonance in a 
two-level system, where the behavior of the susceptibility 
is well understood beyond perturbation theory. Interest- 
ingly, the lowest-order perturbation theory does not ap- 
ply also at exact two-photon resonance, Suj = SuJm-2, as 
discussed below, even though Eq. ([7]) does not diverge. 

III. ANTIRESONANCE OF THE 
MULTIPHOTON RESPONSE 

The field-induced anticrossing of quasienergy levels at 
multiphoton resonance is accompanied by lifting the de- 
generacy of the susceptibilities. It leads to the onset of 
a resonant peak and an antiresonant dip in the suscep- 
tibilities as functions of frequency Suj. The behavior of 
the quasienergy levels and the susceptibilities is seen from 
Fig. [21 For small field amplitude A the multiphoton Rabi 
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FIG. 3: (Color online). Level anticrossing and antiresonance 
of the susceptibilities. The figure refers to a 3-photon reso- 
nance, = 3, in an 5 = 2 system. The involved quasienergy 
states are v = —1 and u = 0. The resonating Zeeman states 
for A = are m = —2 and m — 1 (the ground and 3rd excited 
state). Left and right panels show the quasienergies and 
susceptibilities for the same reduced field A/D. Panels (a) 
and (b), (c) and (d), and (e) and (f) refer to A/D = 0,0.2, 
and 0.3, respectively. 

frequency cx A'^' is small, the quasienergies of inter- 
est £y and £u+i (with 'm{h' -I- 1) — m(i>) = N) come very 
close to each other at resonant 6ui, as do also the suscep- 
tibilities Xi' and Xu+i- 

With increasing A the level splitting rapidly increases 
in a standard way. The behavior of the susceptibilities 
is more complicated. They cross, but sufficiently close 
to resonance they repel each other, forming narrow dips 
(antiresonance) or peaks (resonance). The widths and 
amplitudes of the dips/peaks display a sharp dependence 
on the amplitude and frequency of the field. 

For weak field it is straightforward to find the splitting 
of the susceptibilities 

close to iV-photon resonance between states Ito)*^"^ and 
\m -\- N)^^\ In this region the frequency detuning from 
the resonance 

ALu{m;N) = N{duj - 6ujrn;N)/2 (9) 

is small, |Aci;(m; iV)| < r2fl (m; N). To the lowest order in 
A but for an arbitrary ratio ^/{(m; iV)/|Aa;(m; A^)| the 
quasienergy states and + sgniV) are linear combi- 
nations of the states |m)(°) and |m -I- TV)*^"^. Then from 




A/D 

FIG. 4: (Color online). The multiphoton susceptibility split- 
ting for S = 2. The curves refer to exact A'^-photon reso- 
nances, with A'^ = 1, . . . , 4, for transitions from the ground 
Zeeman state m = —2 to the excited states m — —1, . . . , 2, 
respectively. 

Eq. ([2]) it follows that the splitting of the quasienergies 

AEi^jAT = \ei, - £iy+sgnAf| is 

Ae^.A, = [Acj2(m; N) + 4f7^(m; N)] ^''^ . (10) 

From this expression and Eqs. ([3]), ([5]) it follows that 
the susceptibility splitting is 

Av 8\N\Kim;N) 

A2 [Aio^im; N) + AQl{m; N)]^^^ 

The splitting Axu:N as a function of frequency Slo is 
maximal at A'^-photon resonance, Su — 6LUm;N- The half- 
width of the peak of Axi^-n at half height is determined 
by the Rabi splitting and is equal to V^fln/N. The 
peak is strongly non-Lorentzian, it is sharper than the 
Lorentzian curve with the same half-width. This sharp- 
ness is indeed seen in Fig. [3) Our numerical results show 
that Eq. pT|) well describes the splitting in the whole 
frequency range |Au;| < rtR. 

For small A, the susceptibility splitting is stronger 
than the level repulsion. It follows from Eqs. PTI)) . (fTT|) 
that at exact A^-photon resonance Ae cx A'^' whereas 
Ax cx This scaling is seen in Fig. [4] For 

A/D ^ 1, on the other hand, the eigenstates be- 
come close to the eigenstates of a spin with Hamilto- 
nian —ASx/'2. As a result, the susceptibility splitting 
decreases with increasing A, \Axi^-n\ the propor- 

tionality coefficient here is independent of A^. Therefore, 
for N > 3 Axu displays a maximum as a function of A, 
as seen from Fig. 2] 

A. Two-photon resonance 

As mentioned above, the lowest order perturbation 
theory ([7]) does not describe resonant susceptibility for 
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two-photon resonance. Indeed, it follows from Eq. (fTTj) 
that at exact resonance, Slo = Sio(rn; 2), the susceptibility 
splitting for weak fields is 

Ax.;2 = D-^ [{S ^ m - 1)(S - m) 

X {S + m + l){S + m + 2)]^/^ . (12) 

This splitting is independent of A. The expression for the 
susceptibility ^ is also independent of A, yet it does not 
lead to susceptibility splitting and therefore is incorrect 
at two-photon resonance. 

The inapplicability of the simple perturbation theory 
([7]) is a consequence of quantum interference of transi- 
tions, the effect known in the linear response of multi- 
level systems [iBl . To the leading order in A, the suscep- 
tibility is determined by the squared amplitudes of vir- 
tual transitions to neighboring states. For a two-photon 
resonance, Slo — (5wm:2i the distances between the lev- 
els involved in the transitions Ito)'"' — > |m -t- 1)'"' and 
\m + 2)^°^ |to+1)(°) are equal, e(°)(m-|- 1) - £("'(m) = 
e(°^(m + 1) — e'°)(m + 2). Therefore the transitions res- 
onate and interfere with each other. 

To calculate the susceptibility it is necessary to start 
with a superposition of states |m)^°) and \m + 2)(°\ 
add the appropriately weighted amplitudes of transitions 
|m)(0) -> |to + and |m + 2)^ ^ \m + and 
then square the result. This gives the correct answer. 
The independence of the susceptibility splitting from A 
for two-photon resonance in the range of small A as given 
by Eq. ([12]) is seen in Fig. S] 



IV. SUSCEPTIBILITY CROSSING FOR A 
SEMICLASSICAL SPIN 

The analysis of the simultaneous level and susceptibil- 
ity crossing is particularly interesting and revealing for 
large spins and for multiphoton transitions with large iV. 
For S ^ 1 the spin dynamics can be described in the 
WKB approximation. We will start with the classical 
limit. In this limit it is convenient to use a unit vector s = 
S/S, with s = {sx,Sy,Sz) = (sin cos 0, sin sin (/), cos 6') , 
where 9 and (p are the polar and azimuthal angles of the 
vector s. To the lowest order in equations of motion 
for the spin components can be written as 

Sx = Sy{Sz + fJ.), Sy = -Sx{Sz+ l-i) + fSz, (13) 

sz = -fsy, f = A/2SD, fi = 6lu/SD. 

Here, overdot implies differentiation with respect to di- 
mensionless time r = SDt, that is, s = ds/dr = 
{SD)~^ds/dt. Equations ([13]) preserve the length of the 
vector s and also the reduced Hamiltonian g — H/S^D, 

g = g{9,<P)^-^{s,+tif-fSx. (14) 
For convenience, we added to g the term — /i^/2. 



The effective energy g{6, <j)) is shown in Fig. [5l Also 
shown in this figure are the positions of the stationary 
states s = and examples of the phase trajectories de- 
scribed by Eqs. p^ . 
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FIG. 5: The effective energy g{9, tj)) as a function of the polar 
and azimuthal angles of the classical spin 6 and (p. The lines 
g{9, (f>) = const describe classical spin trajectories. The points 
Ai and A2 are the minima of g, B is the maximum, and S 
is the saddle point. In the region gs > g > gAi there are 
two coexisting types of trajectories. They lie on the opposite 
sides of the surface g{9, (j)) with respect to gs. The plot refers 
to /i = 0.125, / = 0.3. 

An insight into the spin dynamics can be gained by 
noticing that g has the form of the scaled free energy of 
an easy axis ferromagnet [l3|, with s playing the role of 
the magnetization M/M, and with /i and / being the 
reduced components of the magnetic field along the easy 
axis z and the transverse axis x, respectively. In the 
region 

l/p/' + lM^ <1 (15) 

the function g has two minima, Ai and A2, a maximum 
B, and a saddle point S. We will assume that the mini- 
mum A2 is deeper than Ai, that is 

9B> 9s> 9 At > 9A2- (16) 

As seen from Eqs. ([13]) and ([H]) and Fig. [HI for / > the 
minima and the saddle point are located at (j) — and 
the maximum is at (/) = tt; the case / < corresponds 
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to a replacement </> — > + tt. On the boundary of the 
hysteresis region (|15l) the shaUower minimum Ai merges 
with the saddle point S. 

In the case of an easy-axis ferromagnet with free en- 
ergy g, the minima of g correspond to coexisting states of 
magnetization within the hysteresis region (|15|) . For mul- 
tiphoton absorption g is the scaled quasienergy, not free 
energy, and stability is determined dynamically by bal- 
ance between relaxation and high-frequency excitation. 
One can show that, for relevant energy relaxation mech- 
anisms, the system can have coexisting stable stationary 
states inside and outside the region ([TS]). The states cor- 
respond to one or both minima and/or the maximum of 
g; for small damping the actual stable states are slightly 
shifted away from the extrema of g on the (6, (/))-plane. 
We will not discuss relaxation effects in this paper. 

A. Conformal property of classical trajectories 

Dynamical trajectories of a classical spin on the plane 
{6, (j)) are the lines g{9, cj)) — const. They are either closed 
orbits around one of the minima ^i, A2 or the maximum 
B of g, or open orbits along the 4> axis, see Fig. [5j On 
the Bloch sphere = 1, closed orbits correspond to pre- 
cession of the unit vector s around the points Syiij^,s^2, 
or ss , in which s does not make a complete turn around 
the polar axis. Open orbits correspond to spinning of s 
around the polar axis accompanied by oscillations of the 
polar angle 9. Even though the spin has 3 components, 
the spin dynamics is the dynamics with one degree of 
freedom, the orbits on the Bloch sphere do not cross. 

An important feature of the dynamics of a classical 
spin in the hysteresis region is that, for each g in the 
interval {gAi,gs), the spin has two coexisting orbits, 
see Fig. [5j One of them corresponds to spin precession 
around sa^ ■ It can be a closed loop or an open trajec- 
tory around the point Ai on the (0, 0)-plane. The other 
is an open trajectory on the opposite side of the g-surface 
with respect to the saddle point. We will classify them 
as orbits of type I and II, respectively. 

We show in Appendix that classical equations of mo- 
tion can be solved in an explicit form, and the time de- 
pendence s(t) is described by the Jacobi elliptic func- 
tions. The solution has special symmetry. It is related to 
the conformal property of the mapping of Sz onto r. The 
major results of the analysis are the following features 
of the trajectories s(r) of types I and II: for equal g, (i) 
their dimensionless oscillation frequencies 1^(17) are equal 
to each other, and (ii) the period averaged values of the 
component Sx{t) are equal, too, 

<^i{g) = ^iiig), {sx{t))i = {sx{t))ii. (17) 

Here, the subscripts I and II indicate the trajectory type. 
The angular brackets (...) imply period averaging on a 
trajectory with a given g. 

The quantity {sx{t)) gives the classical response of the 
spin to the field oc A. Equation ([TT]) shows that this 



response is equal for the trajectories with equal values of 
the effective Hamiltonian function g. This result holds 
for any field amplitude A, it is by no means limited to 
small A/D where the perturbation theory in A applies. 



B. The WKB picture in the neglect of tunneling 

In the WKB approximation, the values of quasienergy 
£1/ in the neglect of tunneling can be found by quantizing 
classical orbits g{0,(j)) — const, see Ref. [l^ and papers 
cited therein. Such quantization should be done both 
for orbits of type I and type II, and we classify the re- 
sulting states as the states of type I and II, respectively. 
The distance between the states of the same type in en- 
ergy units is fvjj{g)SD p^ . Transitions between states 
of types I and II with the same g are due to tunneling. 

If we disregard tunneling, the quasienergy levels of 
states I and II will cross, for certain values of /z. Re- 
markably, if two levels cross for a given /z, then all levels 
in the range gAi < 9 < gs cross pairwise. This is due 
to the fact that the frequencies uj{g) and thus the inter- 
level distances for the two sets of states are the same, 
see Eq. ([T71) . Such simultaneous degeneracy of multi- 
ple pairs of levels agrees with the result of the low-order 
quantum perturbation theory in A and with numerical 
calculations. 

In the WKB approximation, the expectation value of 
an operator in a quantum state is equal to the period- 
averaged value of the corresponding classical quantity 
along the appropriate classical orbit [Toj. Therefore if 
semiclassical states of type I and II have the same g, 
the expectation values of the operator Sx in these states 
are the same according to Eq. pT|) . Thus, the WKB 
theory predicts that, in the neglect of tunneling, there 
occurs simultaneous crossing of quasienergy levels and 
susceptibilities for all pairs of states with quasienergies 
between gA^ and gs- This is in agreement with the re- 
sult of the perturbation theory in A and with numerical 
calculations. However, we emphasize that the WKB the- 
ory is not limited to small A, and the WKB analysis re- 
veals the symmetry leading to the simultaneous crossing 
of quasienergy levels and the susceptibilities. 

Tunneling between semiclassical states with equal g 
leads to level repulsion and susceptibility antiresonance. 
The level splitting 2flji can be calculated by appropri- 
ately generalizing the standard WKB technique, for ex- 
ample as it was done in the analysis of tunneling between 
quasienergy states of a modulated oscillator [J]- Then 
the resonant susceptibility splitting can be found from 
Eq. ([5]). The corresponding calculation is beyond the 
scope of this paper. 
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V. DEGENERACY LIFTING BY HIGHER 
ORDER TERMS IN & 

The simultaneous crossing of quasienergy levels and 
susceptibilities in the neglect of tunneling is a feature of 
the spin dynamics described by Hamiltonian Q . Higher- 
order terms in Sz lift both this degeneracy and the prop- 
erty that many quasienergy levels are pairwise degenerate 
for the same values of the frequency detuning 5u). The 
effect is seen already if we incorporate the term in the 
anisotropy energy, i.e. for a spin with Hamiltonian 

H^H- lost. (18) 

The Hamiltonian H is written in the rotating wave ap- 
proximation, H is given by Eq. and G is the param- 
eter of quartic anisotropy. The terms S"^ , Sy in the spin 

anisotropy energy do not show up in H even if they are 
present in the spin Hamiltonian Hq but the correspond- 
ing anisotropy parameters are small compared to loq. In 
the rotating frame these terms renormalize the coefficient 
at Si and lead to fast oscillating terms oc S\ exp{±2iujpt) 
that we disregard. 

Multiple pairwise degeneracy occurs where the condi- 
tion on Zeeman quasienergies e^^'^{m) = e^^\m') is si- 
multaneously met for several pairs {m,m'). For G ^ Q 
this happens only for Slo = 0, that is when the modula- 
tion frequency Wi? is equal to the Larmor frequency wq. 
In this case the resonating Zeeman states are |m)'^°) and 
I — m)'-^'' with the same m. The susceptibilities of these 
states are equal by symmetry with respect to reflection 
in the plane {x,y). 

A''-photon resonance for nonzero G and up ^ ujq occurs 
generally only for one pair of states |to)(*'^ and \m+N)'-^\ 
This is seen from panel (a) in Fig. [6l With increasing 
\G\ the difference in the resonant values of frequency in- 
creases, as seen from panel (c) in the same figure. 

The susceptibilities in resonating states are different 
in the weak-field limit. When the frequency ujf adia- 
batically goes through resonance, there occurs an inter- 
change of states, for weak field A: if the state If) was 
close to jm)'^*'^ on one side of resonance, it becomes close 
to \m + Af)*^°^ on the other side. Respectively, the sus- 
ceptibility sharply switches from its value in the state 
|m)'^°) to its value in the state |m -t- iV)'"'. 

Susceptibility switching is seen in panels (b) and (d) 
in Fig. [SI For a weak field the frequency range where 
the switching occurs is narrow and the switching is sharp 
(vertical, in the limit A ^ 0). As the modulation am- 
plitude A increases the range of frequency detuning Slu 
over which the switching occurs broadens. In addition, 
for small G the susceptibility displays spikes. They have 
the same nature as for G — 0. However, they are much 
less pronounced, as seen from the comparison of panel 
(d) in Fig. [H] and panel (f) in Fig. [3] which refer to the 
same value of A/D. 




1,05 
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FIG. 6: (Color online). Quasienergy and susceptibility 
switching for a spin S = 2 with quartic in Sz anisotropy. 
Panels (a), (b), and (d) refer to the dimensionless quartic 
anisotropy parameter G/D = 0.4 in Eq. (|18p . Panels (b) and 
(d) refer to the 3-photon resonance | — 2)'"-' ll)'"' with the 
scaled modulation amplitude A/D — > and A/D = 0.3, re- 
spectively; the dotted line shows the position of the resonance 
5lu/D = 1. Panel (c) shows the dependence of the resonant 
frequency detuning Sujm-.N on the higher-order anisotropy pa- 
rameter G in the limit ^ 0. 



VI. CONCLUSIONS 

In this paper we have considered a large spin with an 
easy axis anisotropy. The spin is in a strong magnetic 
field along the easy axis and is additionally modulated by 
a transverse field with frequency top close to the Larmor 
frequency luq. We have studied the coherent resonant 
response of the spin. It is determined by the expectation 
value of the spin component transverse to the easy axis. 
We are interested in multiphoton resonance where Nujp 
coincides or is very close to the difference of the Zeeman 
energies -E^|.jv ~ i'^ absence of modulation. 

The major results refer to the case where the 
anisotropy energy is of the form —DSl/2. In this case not 
only the quasienergies of the resonating Zeeman states 
|m)(°) and \m + TV)^''^ cross at multiphoton resonance, 
but the susceptibilities in these states also cross, in the 
weak-modulation limit. Such crossing occurs simultane- 
ously for several pairs of Zeeman states. As the modu- 
lation amplitude A increases, the levels are Stark-shifted 
and the susceptibilities are also changed. However, as 
long as the Rabi splitting due to resonant multiphoton 
transitions (tunneling) can be disregarded, for resonant 
frequency the quasienergy levels remain pairwise degen- 
erate and the susceptibilities remain crossing. We show 
that this effect is nonperturbative in A, it is due to the 
special conformal property of the classical spin dynamics. 

Resonant multiphoton transitions lift the degeneracy 
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of quasienergy levels, leading to a standard level anti- 
crossing. In contrast, the susceptibilities as functions of 
frequency cross each other. However, near resonance they 
display spikes. The spikes of the involved susceptibili- 
ties point in the opposite direction, leading to decrease 
(antiresonance) or increase (resonance) of the response. 
They have a profoundly non-Lorentzian shape (jlip , with 
width and height that strongly depend on A. The spikes 
can be observed by adiabatically sweeping the modula- 
tion frequency through a multiphoton resonance. If the 
spin is initially in the ground state, a sequence of such 
sweeps allows one to study the susceptibility in any ex- 
cited state provided the relaxation time is long enough. 

The behavior of the susceptibilities changes if terms of 
higher order in Sz in the anisotropy energy are substan- 
tial. In this case crossing of quasienergy levels is not ac- 
companied by crossing of the susceptibilities in the limit 
A —f 0. Resonant multiphoton transitions lead to step- 
like switching between the branches of the susceptibilities 
of the resonating Zeeman states. Still, the susceptibilities 
display spikes as functions of frequency for a sufficiently 
strong modulating field. 

The results of the paper can be applied also to molec- 
ular magnets in a static magnetic field. The spin Hamil- 
tonian in the rotating wave approximation Q is similar 
to the Hamiltonian of a spin in a comparatively weak 
static field, with the Larmor frequency Suj of the same 
order as the anisotropy parameter D. The susceptibility 
then characterizes the response to the field component 
transverse to the easy axis. Quasienergies £'"^(771) are 
now spin energies in the absence of the transverse field, 
and instead of multiphoton resonance we have resonant 
tunneling. Our results show that a transverse field does 
not change the value of the longitudinal field for which 
the energy levels cross, in the neglect of tunneling. This 
explains the experiment [To| where such behavior was 
observed. 

In conclusion, we have studied multiphoton resonance 
in large-spin systems. We have shown that the coher- 
ent nonlinear response of the spin displays spikes when 
the modulation frequency goes through resonance. The 
spikes have non-Lorentzian shape which strongly depends 
on the modulation amplitude. The results bear on the 
dynamics of molecular magnets in a static magnetic field 
and provide an explanation of the experiment. 

We acknowledge insightful discussions with B.L. Alt- 
shuler and A. Kamenev. This work was supported by 
the NSF through grants No. ITR-0085922 and PHY- 
0555346. 



APPENDIX A: SYMMETRY OF CLASSICAL 
SPIN DYNAMICS: A FEATURE OF THE 
CONFORMAL MAPPING 

Classical equations of motion for the spin components 
(fT5)l can be solved in the explicit form, taking into ac- 
count that = 1 and that g{0,(j)) = const on a classical 



trajectory. For time evolution of the z-component of the 
spin we obtain 

f . r2{ri - rs) - rsiri - r2)sn'^{u;mj) 
Sz(T) = r (Al) 

where ri > r2 > > are the roots of the equation 

[(r + M)2 + 2.g]' + 4/2(r2-l) = (A2) 

and sn(u;r7ij) is the Jacobi elliptic function. The argu- 
ment u and the parameter mj are 

u = LOT, w = ^ [(ri - r3)(r2 - 7-4)]^^^ , 
mj = (ri - r2)(r3 - ri)/{ri - r^){r2 - Vi), (A3) 

Equation (lAip describes an orbit which, for a given g, 
oscillates between = ri and = ^2] the corresponding 
oscillations of Sa;,Sy can be easily found from Eqs. p^ . 

Oscillations of Sz between and for the same g 
are also described by Eq. (|A1[) provided one replaces 
u u + K{mj) + iK'{mj), where K{mj) is the el- 
liptic integral and K'{mj) = K{1 — mj). Clearly, both 
types of oscillations have the same period over r equal 
to 2K{mj)/Cj. They correspond, respectively, to the tra- 
jectories of types II and I in Fig. O that lie on differ- 
ent sides of g{0, (/))-surface. As a consequence, the vibra- 
tion frequencies for the corresponding trajectories ujii{g) 
and uji{g) are the same. This proves the first relation in 

Eq. uni). 

The Jacobi elliptic functions are double periodic, and 
therefore Sz is also double periodic, 

(r) =^ Sz[t + Cj-^{2nK + 2imif')] (A4) 

with integer n, m. Ultimately, this is related to the fact 
that equations of motion (jl3p after simple transforma- 
tions can be put into a form of a Schwartz-Christoffel 
integral that performs conformal mapping of the half- 
plane Im S2 > onto a rectangle on the it-plane. We will 
show now that the mapping has a special property that 
leads to equal period-averaged values of Sx{t) on trajec- 
tories of different types but with the same g. Because 
Sz{t) is double periodic, cf. Eq. (jA4p . so is also the func- 
tion s^{t) = -(2/)-i [2g+{sz{T)+ iif]. Keeping in 
mind that the transformation u u + K{'mj) + iK'{'mj) 
moves us from a trajectory with a given g of type I to a 
trajectory of type II, we can write the difference of the 
period-averaged values of Sx{t) on the two trajectories 
as 

{s^{t))i - {s^{t))ii ^'^^ i s^du (A5) 

where the contour C is a parallelogram on the u-plane 
with vortices at 0, 2A', 3Ar + iK' , K + iK' . It is shown in 

Fig.m 

An important property of the mapping (|XT|) is that 
Sz{t) has one simple pole inside the contour C, as marked 
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FIG. 7: The contour of integration in the u oc r plane. The 
horizontal parts correspond to two trajectories s(r) with the 
same g. The values of s(r) on the tilted parts of the parallelo- 
gram are the same. The plot refers to jj. — 0.125, g — —0.366. 



in Fig. [7l Respectively, Sx{t) has a second-order pole. 
The explicit expression (|Aip allows one to find the corre- 
sponding residue. A somewhat cumbersome calculation 
shows that it is equal to zero. This shows that the period- 
averaged values of Sx on the trajectories with the same g 
coincide, thus proving the second relation in Eq. (|17|) . 
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